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Abstract 

The magnetic phase diagram of a two-dimensional generalized Hubbard 
model proposed for manganites is studied within Hartree-Fock approxima- 
tion. In this model the hopping matrix includes anisotropic diagonal hopping 
matrix elements as well as off-diagonal elements. The antiferromagnetic (AF), 
ferromagnetic (F), canted (C) and paramagnetic (P) states are included in the 
analysis as possible phases. It is found that away from half-filling only the 
canted and F states may exist and AF and P states which are possible for the 
usual Hubbard model do not appear. This is because the F order has already 
developed for on-site repulsion [7 = due to the hopping matrix of the gen- 
eralized model. When applied for manganites the orbital degree is described 
by a pseudospin. Thus our "magnetic" phase diagram obtained physically 
describes how orbital order changes with U and with doping for manganites. 
Part of our results are consistent with other numerical calculations and some 
experiments. 



1 



Typeset using REVT^ 



The three-dimensional (3D) cubic manganites Ri_a;Aa;Mn03 (R is a rare earth ele- 
ment such as La or Nd, A is a divalent alkali such as Sr or Ca) and the layered ones 
La„_„^Sri+„^Mn„03„+i (n = 1, 2) have attracted intensive interest due to their rich mate- 
rial properties Most prominently they show a colossal magnetoresistance which may be 
qualitatively understood on the basis of double exchange (DE) model or the more general 
ferromagnetic Kondo lattice model. However, many recent experimental findings such as 
the doping dependence of charge, orbital, spin ordered phases suggest that additional 
physical mechanisms beyond the DE model should be involved. The possible ingredients 
are: Cg orbital degeneracy electron correlation as well as electron-lattice inter- 

action 0, e.g., by Jahn- Teller coupling P). Though there exist some works to study the 
complex interplay of spin, orbital and lattice degrees of freedom P-[TT[|, it is still very helpful 



to clarify the role of each individual mechanism separately with the assumption that other 
degrees of freedom are frozen out. Actually, experiments showed that the spins in many 
3D manganites Ri_a;Aj.Mn03 are ferromagnetically ordered or A-type antiferromagnetic or- 
dered (ferromagnetic layers coupled antiferromagnetically) in a large region of the doped 
phases, and A-type antiferromagnetic ordered in the undoped phase This means that 



a perfect spin-polarized two-dimensional (2D) plane is always retained over large region of 
doping. Therefore, the spin degree of freedom may be discarded and the following effective 
2D Hamiltonian with only orbital degree of freedom is proposed ||13|,|T^: 



n = -Y. E^^''4c,.' + UY.{n,^ - l/2){n,^ - 1/2) , (1) 

<ij> a a' i 

where < ij > indicates summation over nearest-neighbor sites, t^""' denotes the hopping 
integral and U is the effective inter-orbital Coulomb interaction. The two orbitals — y"^) 
and iSz"^ —r^) have been assumed as pesudospin f and |, respectively. The hopping matrix 
is exphcitly given by tj] = h = 3t/4, t|/ = t2 = t/4, tjj = tj] = = -(+)V3t/4 



along the a;(?/)-direction [|T5|. The most important feature here is the orbital dependence 



of the integral t^?"', which distinguishes the present model from the usual Hubbard model 



where we have = td^a'- 



Although the above 2D model ignores some important effects of three dimensionality 
as we discuss in the later, it may still contain important physics for 3D cubic manganites. 
For this the above model or its strong coupling version, the so-called orbital t-J model has 



been used to study the unusual optical conductivity observed in doped LaMnOs |I^p]6 
At the same time it should be also noted that the above model may be of direct relevance 
to layered manganites, which consist of stacking of single or double Mn02 layers. Thus 
a full understanding of the properties of the above model is necessary. While its optical 
conductivity was studied by others, we discuss a different aspect here. It is natural to ask in 
which properties the present model deviates from the usual Hubbard model. For example, 
as a basic problem, what are the differences between their magnetic (it actually means 
"orbital" for the real manganites) phase diagrams? This will be the main topic addressed in 
this work. Though the orbital ordering for manganites based on the current model has been 



numerically studied at some points [jT3|,|16|, the full phase diagram has not been discussed. 
So the purpose of this paper is twofold: On one hand we want to know how its magnetic 
phase diagram is modified by the generalized hopping matrix, in analogy to the recent 



investigation of the influence of next-nearest-neighbor hopping on the phase diagram. 
On the other hand, as a real model for manganites, our "orbital" phase diagram obtained 
below is certainly useful for the nontrivial question of how the orbital order is changed upon 
doping. 

In the following we use the Hartree-Fock (HF) approximation to study the magnetic 
phases of the Hamiltonian (|l|). Although the HF approximation may overestimate the 
appearance of the ordered states due to neglect of quantum fluctuations, this simple ap- 
proximation may still give many useful information. So it should be kept in mind that the 
quantitative aspects of the results given here should not be taken literally but that these 
results provide a basis for further more elaborated theoretical studies. In this respect the 
so called slave-boson mean field theory is an improved method, but requires much more 
numerical efforts |TB|. Besides the usual ferromagnetic (F), antiferromagnetic (AF) ordered 
and paramagnetic (P) states considered here the canted ordered (C) state is also taken into 



account as a possible phase. Although it has rarely been considered in previous studies of 
the Hubbard model, we include it here motivated by the fact that experiments on LaMnOs 
||T9| have confirmed the type of {Sx"^ — r^)/|3?/^ — r^) orbital ordering which corresponds to 
a canted spin order in the pesudo-spin representation as illustrated in Fig. |1|. Other 
complicated spin ordered states which have been studied for the usual Hubbard model such 



as spiral phase etc. will not be considered. 
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FIGURES 




FIG. 1. Pseudospin representation of staggered |3x^ — r^)/|32/^ — r^) orbital ordering in xy 
plane under the basis: \x'^ — y^) (|) and |3z^ — r^) {[). 

In HF approximation the Hamiltonian (|1]) can be rewritten as (irrelevant constants 
are ignored): 

^HF = - E ll^fc[.c,^' - ^UY^i^S, ■ {S.^ - (S,) ■ (S,)) , (2) 

<ij> aa' i 

where Si = ^ Haa' cla^cra'Cicr' and a are Pauli matrices. We restrict to solutions with uniform 
electron density and canted order of the magnetic moments characterized by 

m(sin6', 0, cos 6') i G A , 
(^.) = <! (3) 
m(— sin6', 0, cos 6*) i G B . 

Here we have assumed a bipartite lattice with sublattices A and B. For later discussion 
it is indispensable to understand the nature of all possible ordered states characterized 
by the parameter 6 (with m describing the size of the order parameter), as well as their 



corresponding orbital ordering |20 



1. ^ = (vr): F order with spin f (J,) or homogenous — y^) (|32;^ — r^)) orbital ordering. 

2. 6 = 71/2: AF order or staggered (|x^ — y^) + |3z^ —r^)) /(|x^ — y^) — |32;^ — r^)) orbital 
ordering. 

3. other 6 values: canted-F or AF order. For example, 6 = tt/S (27r/3) corresponds to 
staggered |3x^ — r^)/|3?/^ — r^) (|?/^ — — z'^)) orbital ordering. 

In addition we note that due to broken SU(2) symmetry by the hopping matrix in Eqs. (1,2) 
the direction of spin polarization in the ordered states is no longer rotationally invariant, 



contrary to the case for the usual Hubbard model. For example, the F phases with spin | 
and I are not equivalent now. We also remark that the spins are lying within the plane in 
the AF ordered state described by Eq. (^. 

With new fermion operators a and h corresponding to sublattice A and B, respec- 
tively, we may write the Hamiltonian in the momentum space as T^hf = 'iNUm? /?> + 
('^It ^It '^fci ^Ii)'^^('^fcT ^fc? '^fci ^ki)^ with T-C^ denoting a 4 x 4 Hermitian matrix. Here 

is the total number of lattice sites and the reduced Brillouin zone is l/c^^l + l^j^l < tt or 
— 7r/2 < fci, k2 < 7r/2 {k^^y = ±ki + ^2). In principle leads to four energy bands 
(a = 1,2,3,4). The ground state energy is given by Eqs = J2k,a^k + 2A^f/m^/3, where 
the first contribution depends on the band filling n = N^/N where A'e is the total electron 
number. Next we calculate numerically the values m, 6 by minimizing the Eqs for given 
repulsion U and band filling n or doping 6 (= 1 — n). For this purpose we have used a mesh 
of 10^ X 10^ k points in the Brillouin zone. In the following discussion of results t is set as 
energy unit. 

With values of m, 6 obtained from the minimization of Eqs we may plot the magnetic 
phase diagrams in the plane {6, U) , as shown by Figs. and |^. Although there are a lot of 
studies for the usual Hubbard model (i.e., ^1 = ^2 = 1, ^3 = 0) we give its phase diagram 
again in Fig. ^ since the canted ordered state is added here. We also want to compare it 
with the result for the current generalized Hubbard model whose phase diagram is displayed 
in Fig. ^ (solid line). In the present model, the hopping matrix involves anisotropic diagonal 
elements and off-diagonal ones. To see more clearly how the phase diagram is influenced by 
them, we plot it for a simplified model with a hopping matrix which only includes anisotropic 
diagonal elements (i.e., ti = 3/4, t2 = 1/4, ta = 0); it is denoted by a dashed line in Fig. |^. 



6 




0.1 0.2 0.3 0.4 0.5 

Doping 6 

FIG. 2. Magnetic phase diagram for the usual Hubbard model with inclusion of ferromagnetic 
(F), antiferromagnetic (AF), canted (C) ordered and paramagnetic (P) states. 




0.1 



0.3 



0.4 



0.2 

Doping 6 

FIG. 3. Magnetic phase diagrams. The solid (dashed) line are for the generalized Hubbard 
model with h = 3/4,^2 = 1/4, ta = \/3/4 {h = 3/4,^2 = 1/4, ts = 0) respectively. F(P) denote the 
stable phase at \J=6=0 for both cases. 



/^From Figs. ^ and 0, a number of differences can be seen. Firstly, the phases for the 
usual Hubbard model include AF, F, P states as well as a narrow region with a canted 
phase. It appears in a crossover region from AF to F phase, which is understandable. The 
phase diagram is then qualitatively modified by introduction of anisotropic diagonal elements 
where only the canted and F ordered phases exist and the original AF and P phases are 
excluded away from half- filling (i.e., 6 > 0). Furthermore when the off-diagonal elements are 
introduced the phase diagram is mainly quantitatively modified, i.e., the stability region of 
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the F ordered phase is enlarged. We point out that the F phase in Fig. ^ always means that 
spins align in the up direction (|) due to a larger hopping matrix element t]j (i.e., ti > t2). 

The above results may be understood if we consider the property of the kinetic term 
alone, i.e., assume that U = 0. Then the Hamiltonian (|I]) may be exactly diagonalized with 
two bands: 

= — (ti + t2)(cos kx + cos ky) ± ^ {ti — t2)^(cos kx + cos ky)"^ + 4it|(— cos + cos ky)"^ . (4) 

^From the 2;-axis magnetization = (c,^|Ci| — c]^|Qj)/2 may be calculated, as shown in 
Fig. 1^ where the curves Mz/n vs 6 are given. For the usual Hubbard model Mz is always 
zero, corresponding to P state for the abscissa in Fig. However, a finite magnetization 
(with spin |) is always exhibited for the present model as shown by the solid line in Fig. ^; 
also for the simplified model (zero off-diagonal hopping) shown by the dashed line, except at 
the point 6 = 0. Moreover, the Mz value for the current model is always bigger than that for 
the simplified one, which means that the F order of the former is stronger. Because of the 
underlying F order already in the U = case only the canted (with nonzero F component 
involved) and F ordered phases may exist in finite U region away from half-filling. It is 
also easily understood why the F phase region is enlarged after inclusion of the off-diagonal 
hopping matrix elements. 

0.175 
0.15 
0.125 
^ 0.1 
S 0.075 
0.05 
0.025 

0.1 0.2 0.3 0.4 

Doping 6 

FIG. 4. The z-axis magnetization for the generalized Hubbard model (solid line: = \/3/4, 
dashed line: ts = 0) for on-site repulsion U = 0. 
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We consider the case of half-filling (i.e., 6 = 0). It has been shown that the usual 
Hubbard model has an AF ground state for any finite U, corresponding to the ordinate 
in the phase diagram of Fig. 0. We speculate that this conclusion might also be true for 
our present generalized Hubbard model. Although the present ground state is F ordered at 
t/ = as given in Fig. ^, it is found to have AF order for U > 1.2. This was concluded 
from a calculation employing a mesh of 10^ x 10^ k points. Whether or not there exists a 
crossover between F and AF for finite U is not yet certain. 

Our "magnetic" phases obtained from Hamiltonian ([l|) actually correspond to the orbital 
orderings relevant for the manganites as mentioned before, therefore we now discuss them 
in detail. To see clearly how the orbital order changes with doping 6 we have plotted in Fig. 
^ the relation of parameter 6 with 6 for fixed U. It is found 6 decreases monotonically from 
tt/2 to with increase of 6 indicating the process AF— C— F order for the pseudospins 
as shown in Fig. ^. Let us take ?7 = 10 as an example in the following. (The real value 
for it in manganites is about 4 ~ 10 In the undoped case, 9 = 11/2 means a staggered 
(|x^ — y"^) + |32^ — r^))/(|x^ — y^) — {Sz"^ — r^)) orbital ordering, which is consistent with 
the result by Motome and Imada from quantum Monte Carlo calculation [113, but slightly 



different from the experimentally observed |3x^ — r^)/|3?/^ — r^) This may be due to 
some other physical mechanism, e.g, a Jahn- Teller effect which is present in real materials 
but not included here [|l^]. On the other hand, it is found that the |3a;^ — r^) /|3?/^ — r^) type 



orbital ordering may be also exhibited from our pure orbital model at very small doping 
6 ~ 0.06. With increase of doping the orbital order changes continuously until at some 
critical value Sc — 0.12 which is still small it becomes — y"^) type. This behaviour for 
small 6c has also been noted by Horsch et al. based on the orbital t-J model (i.e., large U 



version of ours) with the exact diagonalization method [|T^]. Interestingly it is also seen that 
the value 5c is not a monotonic function of f/, which reaches a maximum ~ 0.2 at about 
U = 5. Recently Akimoto et al. found an A-type antiferromagnetic metallic ground state 
for La1_j.Srj.MnO3 at doping x > 0.5, and they proposed that — y"^) orbitals should be 
occupied in this state This is consistent with our result derived from the planar model 
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(0). In addition, we also noticed the experimental findings that many cubic manganites are 
3D ferromagnets in the intermediate doping regime nearby 6 = 0.3 [p!2|,p|, where the hopping 
along the 2;-axis is allowed by the DE mechanism. Thus it is instructive to discuss the effect 
of three dimensionality neglected in our model in this doping region, see also Ref. [0. For 
a 3D model, the hopping along the 2;-axis is orbitally highly anisotropic; only the hopping 
between two |3z^ — r^) orbitals is permitted. So the 2D ferro-orbital order with — y"^) 



type polarization, as found above, may be suppressed in order to gain the kinetic energy 
along the 2;-axis. Such a depolarization may be simulated by introducing a negative chemical 
potential /i on the |3z^ — r^) orbital ||14[. It is found that, with reasonable fi valid for 3D 
manganites the disordered orbital liquid state may appear at some doping nearby 6 = 0.3. 



This orbital liquid was first discussed by Ishihara et al. , which may be enhanced by the 



quantum fluctuations beyond HF approximation p5[. On the other hand, for the layered 
manganites with single or double layers of Mn02 where the inter-layer or bilayer hopping (or 
correspondingly /i) is assumed very small, the orbital liquid state may be excluded and our 
ferro-orbital ordered state is still expected to exist in the doping region mentioned above. 
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FIG. 5. Tlie variation of parameter 9 (in units of vr) characterizing the orbital ordering with 
doping 6 for fixed U for the generalized Hubbard model (|l]). 

In the above discussion we have shown the possible orbital ordering as function of doping 
5 and repulsion U . The results are partly consistent with previous numerical calculations, 
and also with some experimental results on manganites although only orbital degree of free- 
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dom is retained in the present model. For a full understanding of various phases observed 
in manganites experimentally, one needs to consider the spin and orbital orderings simul- 
taneously [jlO|,|lT|, as well as include lattice degree of freedom in some cases |T^. However, 
our orbitally ordered structures obtained above are very instructive to understand the role 
of orbital degree of freedom itself, i.e., how the interplay of the orbital kinetics (with atten- 
tion to its peculiar matrix character) and orbital correlation induces the orbital ordering. 
Moreover, they may be also possibly observed under some conditions, e.g., at large external 
magnetic fields. 

In summary, within HF approximation we have studied the magnetic phase diagram of 
a generalized Hubbard model proposed for manganites where the hopping matrix includes 
anisotropic diagonal elements and off-diagonal elements. The AF, F, canted ordered and 
P states are considered as possible phases. It is found that in the doped case only the 
canted and F states may exist, while the AF and P phases of the usual Hubbard model are 
not favored. This is because the F order has already been established at f/ = due to the 
peculiar hopping matrix. When applied for manganites, the spin index of this model actually 
represents the orbital degree of freedom. Thus our obtained magnetic phase diagram is able 
to describe how the orbital order changes upon doping. 

Two of authors (Q. Yuan and T. Yamamoto) would like to thank the support of Visitor 
Program of MPI-PKS, Dresden. Q. Yuan also acknowledges the part support by Chinese 
NSF. 

Note Added in Proof: Recently the scenario for phase separation has been extensively 
studied for various models in the manganites ||26|. Such a tendency of separation into 



orbitally ferro- and antiferro- ordered regions may also exist in our case of the canted phase 
in Fig. 3. This problem deserves further detailed investigation. 
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